General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FACILITY FORM 602 


77H S' <9 7 6-0 


(} 


i — ; — v> 


A COMPARISON OF HUMAN RESPONSE MODELING 
IN THE TIME AND FREQUENCY DOMAINS 


By Lawrence W. Taylor, Jr, 
NASA Flight Research Center 
Edwards, Calif, 


GPO PRICE $„ 
CFSTI PRICE(S) $_ 


Hard copy (HC) 
Microfiche (MF) 


IlM 

/ As"' 







Presented at the 

USC/NASA Conference on Manual Control 
Los Angeles, Calif. 

March 1-3, 1967 


FkepG/AJT 


A COMPARISON OF HUMAN RESPONSE MODELING 

IN THE TIME AND FREQUENCY DOMAINS 

By Lawrence W. Taylor, Jr. 

NASA Flight Research Center 
Edwards, Calif. 

SUMMARY 


Frequency and time domain methods of analyzing human control response while 
performing compensatory tracking tasks are reviewed. Sample linear model results 
using these methods are compared and discussed. The inherent requirement of con- 
straining the freedom of the form of the pilot models is also discussed. The constraint 
in the frequency domain consists of smoothing with respect to frequency; whereas, the 
constraint for the time domain model is more natural and meaningful in that it consists 
simply of limiting the memory of the pilot model. The linear models determined by 
both methods were almost identical. 

The time domain method of analysis enables the determination of a nonlinear pilot 
model. The inclusion of a cubic as well as a linear term accounted for only a small 
additional part of the pilot’s remnant and indicated that only a small portion of the total 
power of the pilot's output is caused by nonlinearities. The power-spectral density of 
an ensemble average of the pilot's output is used to determine the upper limit of the 
amount of power associated with a deterministic response. The indication is that more 
than half the remnant is stochastic when a linear model is used. 


INTRODUCTION 


Since the MIT-NASA Working Conference on Manual Control at Cambridge, Mass. , 
February 28-Mareh 2, 1966, additional work in human response analysis, both theoreti- 
cal and experimental, has been performed at the NASA Flight Research Center. Much 
of this work has been part of a continuation of a NASA-USAF-Cornell Aeronautical 
Laboratory program to obtain pilot describing functions from both flight and simulator 
tests. In addition to this program, the Flight Research Center has initiated a study, 
under the guidance of Dr. A. V. Balakrishnan of UCLA, of nonlinear time domain 
methods as applied to the problem of modeling the pilot in a compensatory tracking task. 

It is the purpose of this paper to assess first the frequency domain method of 
analysis and then the time domain analysis. A comparison of the results of the two 
forms of analysis applied to a linear model is made, and their advantages and disadvan- 
tages are discussed. Next, the time domain method of analysis is applied to the identi- 
fication of a nonlinear pilot model. This is the first time that the nonlinear time domain 
method has been applied to human response data. The power-spectral density of an 
ensemble average of the pilot's output is used to estimate the amounts of power that 
correspond to linear, deterministic, and stochastic control response. 



DESCRIPTION OF EXPERIMENT 


The classical experiment for obtaining data from which pilot models can be identi- 
fied is illustrated in figure 1. The pilot is asked to minimize the error, e, displayed 
to him by an oscilloscope, television screen, or meter by manipulating a controller. 

The controller deflection, c, is sent to an analog computer which computes the response 
of the controlled element and adds to it the input disturbance function, i, forming an 
error which, in turn, is sent to the display. The signals are either processed during 
the experiment or recordings are made of the signals which are later processed to obtain 
the model of the pilot (ref. 1). Similar experiments have been performed in flight in t 

which the pilot maneuvers the airplane (refs. 2 and 3). Most of the data analyzed in this 
paper were collected as part of the joint NASA-USAF-Cornell Aeronautical Laboratory 
human response studies (ref, 3) involving the T-33 variable-stability airplane and ground » 

simulators , 


DISCUSSION AND RESULTS 


Frequency Domain Methods 

Classically, the model of the pilot is considered to be a linear-describing function 
with output, o, plus a remnant signal, r, as shown in figure 2. The describing func- 
tion, Yp(jw), can be obtained by first computing the cross -spectral density functions 

§i c (ju>) and $ie(j w )* The estimate of Yp(jcj) is then given by the ratio (from ref. 4) 


Yp(jca) = 




Cross -spectral density functions have generally been used instead of Fourier 
transforms (ref. 5) as a means of removing the bias in the estimate of Yp(joj) intro- 
duced by the remnant. The use of cross-spectral density functions, however, was 
shown in reference 6 to have no effect on the bias. The same estimate of Yp(ju>), 
therefore, can be expressed as the ratio of Fourier transforms 


and the bias in both cases is 



F[c(t)l 

F[e(t)] 


Vi") - Y P<i“> = 


$ir(i M ) 


Ffr(t)1 

F[e(t)] 


These conclusions, which were developed in reference 6, depend on the ability to 
express the cross-spectral density in terms of Fourier transforms. Appendix A con- 
tains a detailed mathematical proof of the expression, and appendix B presents experi- 
mental results which further substantiate these conclusions. 
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The frequency domain method of analysis was used in the human control response 
study of reference 3 in which compensatory tracking experiments were performed in 
flight and on simulators. Figure 3 shows some simulator results in which the values 
presented for 1?p(jw) are the means of 10 runs for pilot A and three runs each for 
pilots B and C. The vertical lines indicate the range of plus or minus one sigma for 
each of the points. The lack of a vertical line indicates the range to be less than the 
height of the symbol. The input disturbance function consisted of the sum of 10 sinusoids. 
Values of 1fp(ju>) were determined at the input frequencies. The use of sinusoids for 
the input disturbance function has the advantage of concentrating the power at several 
frequencies, thereby enhancing the accuracy of the estimate of the pilot describing func- 
tion at the input frequencies . 

If, on the other hand, a random input is used, mathematical difficulties may be 
encoimtered when using frequency domain methods. For example, if no constraint is 
placed on the form of Yp(jw), the resulting estimate will account for the entire pilot 
output, c, erroneously indicating the remnant, r, to be zero. This result comes from 
applying the relationship 


0 


D 


tp(jw) = 


§ie« w ) = ne<t)] 


at all frequencies. Figure 4 shows an example of the erratic function (dashed lines) 
that results from such a procedure. A constraint on Yp(jw) is usually provided by 
smoothing the values of the cross-spectral density functions or Fourier transforms, or 
by fairing a curve through the raw estimates of the pilot’s describing function, 'Yp(jw), 
or both. Nevertheless, the required constraint on “Yp(jw) compromises one of the 
claims made for the analyses in the frequency domain, namely, that of unlimited free- 
dom in model representation. Also indicated in figure 4 is the effect of smoothing on 
the average linear coherence, p a . The value is seen to change from 1. 0, for the raw 

estimates, to 0. 91 for the faired case. The value of p a or any other measure of the 
remnant has meaning only when connected to some particular pilot-constrained model. 

If the raw cross spectral densities are smoothed, the estimate of Yp becomes 


= 


X Wk$i C (jw + jkAcu) 
k=-K 


k 2 K W k$ie( ],cu + 31 ^ Aco ) 


This estimate can also be expressed in terms of Fourier transforms as 

X W k F k *[i(t)] F k [c(t)] 

^ 

£ W k F k *[i(t)]F k [e(t)] 
k=-K 
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It should be noted that although the smoothing is indicated to be in the frequency 
domain the same result can be obtained by averaging or truncation in the time domain. 

Another point made in reference 6 is that the correlation coefficient, as it is 
usually defined, always has the value of unity 



This result is now recognized to be true in the absence of smoothing. It is now sug- 
gested that the power-spectral density be smoothed over a finite frequency bandwidth. 
The modified definition of p then takes the form 


p 2 >) = < 


2K + 


t 2 


k~K 


+ jkAw) 


2K + 1 


.§ii("+ ^ w >2frrT,A $cc( w + 


k=-K “* T ‘S 

V[i(t)]F kt o(t)]| 

ks-K 


f F k * [i(t)] F k [i(t)] f F k *[c(t)] F k [c(t)] 
fc=-K k£-K 


^ 1 


The linear correlation coefficient so defined is not, in general, equal to unity even for a 
1 inear system, contrary to popular belief. The coefficient does equal unity if the rela- 
tionship between i and c is constant over the bandwidth, whether or not the system is 
linear. It is apparent, therefore, that a further improvement is needed. It is suggested 
that the linear correlation coefficient be expressed as 


Q 


ft 


p 3 2 (co) = 


Vr*[i n (t)](F[o n (t)]) 


F*[in<t)]F[i(t)] X F*[c n <t)]F[c n (t)] 


n=i 


In this expression, an ensemble average is used instead of an average over a range in 
frequency. The result is that now only linear systems free of any stochastic signals 
will produce a value of unity, provided the input is not the same for each ensemble. In 
the event the input is the same, nonlinear systems will appear to be linear, with p 3 = 1. 


ft 


» 


Linear Time Domain Method 


Let us now consider a linear analysis in the time domain in which the output of a 
linear pilot model is expressed in the form (see ref. 7) 


r T M 

c(t) = J hp(r)e(t - r)dr 
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Because the time histories c(t) and e(t) must be sampled for analysis, it is more 
appropriate to write 


c(n) = 



hp(m)e(n - m + 1) 


or in matrix form 


c = Ehp 


where 


'e(M) e(M - 1) , , . e(3) e(2) e(l)’ 
. e(4) e(3) o(2) 

e(4) e(3) 
e(4) 

e(N - 1) e(N - M) 

,e(N) e(N - 1) . . . e(N - M + 1). 



The sampled impulse response of the pilot model, hp(m), can be obtained by using 
least squares 

h p = [eTeHeTc 

Inherent in the time domain representation of the pilot model is the assumption that 
the output at any one time is a function of only a finite time of the history of the error. 
The finite time period (or maximum memory) is denoted by T]y[ (T^ = MAr) in the 

integral expression or M in the summation expression for the pilot model output. For 
the pilot model, T]y[ was varied (by changing At and keeping M constant) until it 
was determined that the value of hp(r) was essentially zero beyond about 1 second. 

The value of Tjyj selected is somewhat larger than 1 second. Figure 5 shows an 
example result of such an analysis. It can be seen that the model impulse response 
first peaks at about 0.25 second, then reverses at about 0. 45 second to peak in the 
opposite direction at about 0. 6 second, then subsides to zero. The first sample 
(r = 0. 05 sec) is typically negative but has been faired to correspond to a pure time 
delay of 0. 05 second. One indication of the degree to which a model represents an 
actual pilot is the ratio of the output of the power of the model in relation to the total 
power of the pilot's output. Linear pilot models will typically account for 65 to 90 per- 
cent of the total power of the output for a 4-minute run. The percentage is somewhat 
higher for shorter runs. 
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The time domain results can be transformed to the frequency domain for compari- 
son with the frequency domain results through the use of the Fourier transform 

$p(jW) = p[fip(T)] 

Figure 6 shows such a comparison. The agreement between the two methods is good, 
One advantage of the time domain method is that Xp(jw) can be determined as a con- 
tinuous function of frequency even when the input consists of sinusoids as it does for the 
example shown. The time domain results by their nature will always yield ^Yp(jco) = 0, 

o d lVH 

or -180 , as oj — 0 and — L =0, as w— ► 0. Had the pilot describing function 

exhibited gain fluctuations or a phase lag at the lowest frequencies shown, these charac- 
teristics would not have been identified by the time domain model unless the value of 
T M were increased considerably. This does not mean that the time domain model 

causes a restriction, but, rather, that T^j should not be unduly limited. Just as it is 

necessary to constrain the frequency domain model by smoothing, it is also necessary 
to constrain the time domain model by limiting Tm to a value considerably less than 

the record length. For pilot models, this represents a very natural and meaningful 
constraint in the time domain, compared with smoothing in the frequency domain. Still 
another advantage of analysis in the time domain is the capability of constructing non- 
linear pilot models. 


Nonlinear Time Domain Method 

Nonlinear behavior on the part of the pilot accounts for at least part of the remnant 
of a linear pilot model. It is, therefore, of interest to investigate nonlinear pilot 
models. The output of a nonlinear time domain pilot model can be expressed by using 
a Volterra integral series 


r T M 

c(t) = J h Pl (r)e(t - T)dr 


(linear) 


+ 



h P2 (T l’ T 2 )e(t " T l )e(t " r 2> dr l dr 2 
(quadratic) 



» r 2’ r 3> e ( t - T ])e(t - r 2 )e(t - 73^7-^7^73 
(cubic) 


+ •«»#••••«• 
(higher order) 
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or in the discrete case 


c(n) * ^ h|(m)e(n - rn + I) 
ni=l (linear) 

M M, 

+ y V h2(mi> m 2 )e(n - mj. + l)e(n - m 2 + 1) 
m l”l m 2 = * (quadratic) 

^ M M 

+ ^ ^ 2^ hg(mi, m 2 , 1113)0(11 - ni| + l)e(n - m2 + l)e(n - m3 + 1) 

m^l m 2 -l m 3 =l (cubic) 

+ »#*»»»#*»»♦ 1 t 

(higher order) 


It was reasoned that the pilot's control response would be symmetrical so that only the 
first (linear) and third (cubic) terms were used. The algorithm used to perform the 
analysis was again based on a least-squares solution (ref. 7). 

If 


where 


E 1>s h 



c 




/hid) \ 

/ hi<2) 


h x (M) 

h 3 (l. 1, 1) 

had. 1,2) 

* 

1*3(1, 1 , M) 

1*3(1, 2, 2) 


(hg(M,M,M), 
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e(M) . . . , e(l) (e(M)e(M)e(M)) . . (e(M)e(M)e(l)) (e(M)e(M - l)c(M - l)) . 
e(M+ 1) . . c(2) (e(M)e(M - l)e(l)). .(e(l)o(l)o(l)) 

c(M 4- 2) 
e(M + ») 


, (e(N - M + l) 2 e(N - M)) 

e(N) . . . , o(N - M + 1) (e(N)o(N)e(N - 1)) e(N -M+ 1)3 


Then 



It is difficult to present the results of such an analysis in a meaningful form, ?..it it 
is instructive to look at an example step response. Figure 7 shows the response to (1) a 
step of very small amplitude so that only the linear term contributes significantly to the 
response, and (2) a large step. The responses have been normalized to the amplitude 
of the step inputs to facilitate comparison. The response of the pilot model to the 
larger step is slightly faster, has more overshoot, and has a lower steady-state gain. 
The inclusion of the cubic term increased the ratio of the power of the model output to 
the total power of the pilot output by only a few percent. This result would indicate the 
remnant to be largely stochastic as opposed to nonlinear and deterministic. 

If the nonlinear model is expanded to include more samples of the cubic term and 
higher order terms, dimensionality will become a problem. One means of reducing the 
total dimension is offered by Balakrishnan (ref. 7) and Hsieh (ref. 8). Through the 
employment of the adjoint system of equations, a cubic weighting function of the form 

L 

%(t - r)f 3 (t - T 2 )%(t - T3)dr 1 dr 2 dT3 


is obtained. This equation reduces the weighting function of three variables to a single 
function of one variable, This technique has not yet been applied to the problem of 
identifying pilot models, and it is not known if the reduction in dimensionality justifies 
the added conputation required. 


Analysis of the Pilot’ s Output 

It would be of interest to know what portion of the pilot’s response is deterministic, 
but not linear, in order to assess the potential of a nonlinear pilot model in describing 
the pilot’s output. It is blown that at least part of the pilot’s output is stochastic, since 
results of repeated experiments are never identical. It is possible to estimate the pro- 
portioning of the power of the pilot's output by examining the power -spectral density 
functions of both the pilot’s output and its ensemble average. Both functions are shown 
in figure 8, The cross-hatched peaks in the graph show the amount of power associated 
with a linear response at the frequencies of the input. The shaded areas show the change 
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in the power as a result of ensemble averaging. Since the deterministic response would 
be unchanged by averaging, the shaded areas are an indication of the power associated 
with the stochastic portion of the pilot's output, which will not be accounted for by a 
deterministic model. The unshaded areas, then, are upper limits on the potential 
increase in power accounted for by using a nonlinear rather than linear pilot model, 

The bar graph at the right of figure 8 shows the proportioning of the power of the 
pilot's output to 1)0 91,7 percent linear (and time invariant), 4. 5 percent stochastic, and 
3, 8 percent nonlinear and other types of responses. These results should not be gener- 
alized, since changes in the controlled element and input can cause a marked change in 
the proportioning of the power. It should also be noted that a small amount of power 
may be attributed to a nonlinearity that is significant in other aspects, such as limit 
cycles , 


CONCLUDING IDS MARKS 


A review of frequency and time domain methods of analysis shows that both methods 
require constraints on the freedom of the pilot models. The constraint in the time 
domain is more natural and straightforward than that of smoothing in the frequency 
domain. The two methods show good agreement for the linear model when the input 
disturbance function consists of sinusoids. 

The inclusion of a cubic term in the time domain pilot model represents the first 
time the analysis has been applied to human response data. For the example discussed, 
only a few additional percent of the power of the remnant was accounted for by the addi- 
tion of the cubic term. An investigation of the power-spectral density of an ensemble 
average of pilot output indicates the reason to be the largely stochastic nature of the 
remnant. The proportioning of the power of the pilot's output appears to be about 
92 percent due to linear response, 4 percent due to stochastic response, and 4 percent 
due to nonlineur and other types of responses . 

With this step toward the application of time domain methods of analyzing human 
control response, the work ahead holds much promise for the determination of more 
meaningful and useful pilot models. 
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APPENDIX A 


COMPARISON OF TWO ESTIMATES OF CROSS-SPECTRAL DENSITY 

APPLIED TO RANDOM SIGNALS 


Two samples of a random signal were used to test the equivalence of two estimates 
of cross-spectral density. Figure 9 shows time histories of the two samples for 250 of 
the 400 time points used. The cross -correlation function given by 


%y(r) 


L 

2T 


/ 

-T 


x(t)y(t + r)dt 



x nyn+m 


where 


r » mh 


h = time interval 


was computed and is plotted in figure 10. The random nature of the signals is borne out 
by the erratic nature of the cross-correlation function. 


Figures 11(a) and 11(b) show the real and imaginary parts of the cross-spectral 
density as estimated by two different expressions, as follows: 


$xy^0 w ) ” * pxy( T )] 


h ^ 
N+l 

k=^N 



x nyn+k 


and 


txy 2 0“) = 


F*[x]F[y] 

2T 






The values obtained by using the preceding expressions were identical except for an 
occasional difference in the fourth significant figure, 


One possible explanation why calculations made by other investigators resulted in 
a difference for the two estimates is the effect of using a limited number of values (lags) 
of the cross -correlation function. Figures 12(a) and 12(b) show the effect of using fewer 
lags than the maximum number possible (in this case, 400). As fewer lags are used, 
the estimate departs from the example for 400 lags, nullifying the equivalence of the 
two estimates. 


Another possible explanation for an apparent difference in the two estimates could 
be the effect of smoothing shown in figures 13(a) and 13(b) , If smoothing were used for 

$xy^(i 0J ) = F[R xy (r)] and not for $ xyo (jw) = , the difference could 
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APPENDIX A 


erroneously be Interpreted as being caused by the form of the estimates of cross- 
spectral density. 

In conclusion, the two forms of cross -spectral density estimates yield identical 
results if (1) the number of "lags” is not limited, and (2) if identical smoothing is used 
in both bases. 
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API 5 ENDIX B 


EQUIVALENCE OF TWO ESTIMATES OF CROSS-SPECTRAL DENSITY 
Consider the following sample time histories: 



x(t) = 0 for t > T and t < -T 
y(t) = 0 for t > T and t < -T 

An estimate of the cross-correlation function can be expressed as 

T 

+ R xy (r) = jy f x(t)y(t + r)dt 
and the Fourier transforms by 

T 

F[x] = f e"'3 aJt x(t)dt 
-T 

T 

F[y] = f e"3 cut y(t)dt 
-T 


4 - 

Another estimate of R X y(r) involves division by 2T - r instead of 2 T, but the 
expression shown is usually more desirable. 
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APPENDIX B 


The most popular estimate of the cross-spectral density function is the Fourier 
transform of the cross-correlation function 

2T T 

$x yi a«) = FpVrt] = .4 e ' iWT W_^ *<t)y(t + T)dtdT 

An alternate and equivalent estimate is the product of two Fourier transforms 

$xy 2 <3 w ) = = W Jp x ( t ) e: ’ a,tdt y(«)e" 3C0Q! dQ! 

in which F*t ] denotes the complex conjugate of Ft ]. 

It is now desired to show that the two estimates are equal. This is done by sub- 
stituting a change in variables and rearranging terms, Let a = r + t and substitute 
r = a - t in the expression of the first estimate 

2T+t T 

fcy^) = £t/ / e"3 aj ( a “ t )x(t)y(oi)dtdQ! 

— 2T+t — T 

Examination of the values for which x(t) and y(r) and y(a') are zero enables the 
limits of integration to be simplified as illustrated in the following figure. Before the 
change of variable, the combinations of t and f for which the integrand is not 
necessarily zero are represented by the shaded area. Because x(t) and y ( t) are zero 
for t, r > T, or > -T, the integrand is necessarily zero outside the shaded area. 


Original limits of integration 


Integrand is zero outside this area 
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After the change in variable, the same area is transformed to that shown in the following 
figure : 



Superimposed on the same plot are the simplified limits of integration. These limits 
have no effect on the integral, since both sets of limits completely cover the area for 
which the integrand is not zero. Therefore 

T T 

$ X y 1 (j«)=^r/ f e-Ma-^xWy^dtda 
It is now possible to rearrange the expression to get 




-T . 

/ e"^ ojQi y(a)dQ; 

-T 


This expression is recognized to be the second (postulated) estimate. Thus 

$xy^ = 2T F *[ x l F [y] = $xy 2 
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appendix b 

E cross- (or power) spectral density are, therefore, equal. 


The two estimates of 

speoteal'denslties Is offered'for the" sampled time histories where 


A similar proof of the equivalence of : the two > of cros S . end power- 


r»T 

-T 

F[y) = f y(a)e-i w “d“ » m jfl /2 yme-i<umh 

The second estimate of the oross-speetral density is then 


jeonh 


F*M = l x(t)ei“*dt « 


$ X y 2 (j" ) =^pF*M F [y] K iS4 


.n=% 


x n e 


j co nil 


LniA/2 


y m e 


-jeomh 


The first estimate may be written in the form 


$ xyi (jo>) = F[R x y(T>] - k 2^ e i Idl N+l n jf^/2 


x n y n +k 


Then, letting m = k + n, and substituting k m n 


N+n 


$xy 1 ^ a) ^ = N+l j- 


e 


-jco(m-n)h 


$4? 


m=-N+n 


n=ls/2 Xn ' 


y- 


m 


Rearranging 

X ,■ ,_j!_ ^ e)“nh x ^ e -i" mh y m 
$xyi<l w ) N+1 n =tft/2 “m^J+n 

1 -m n > n/ 2 and < -N/2, the limits can be changed 

Since x n and y m equal zero for m, n^JN/z, an 

to ±N/2. Therefore 

ft / 1w \ Ik- e3 Wnh x n ^ e -jcomhy _ ^ 

$x yi <3") " N+1^4/2 n m =^/2 
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SYMBOLS 

pilot output (control deflection), inches 

mean square or total power of c, inches 2 

error matrix 

error, radians 

Fourier transform 

time interval, seconds 

impulse response of pilot, inches/radian 

input (external disturbance function), radians 


maximum value of k 

index for frequency 

total record length, seconds 

maximum value of m, M = 

Index for the argument of h p 
maximum value of n 
index for time 

linear output of pilot model (control deflection), inches 
cross -correlation function 

remnant signal of pilot model (control deflection), inches 
Laplace variable 

one-half total record length, seconds 

maximum memory time of the pilot model, seconds 

time, seconds 
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H 


W k 

x»y 

Y c (jw) 

Yp(jco) 

a 

P 

Pa 

r 

At 

<M w ) 

$xx(“) 

00 

Aco 


weighting function 
sample signals 

controlled element transfer function, radians/inch 

pilot describing function, inches/radian 

time variable, seconds 

linear correlation coefficient 

average linear coherence 

argument of hp, seconds 
incremental value of r, seconds 
cross -spectral density of x(t) and y(t) 

power-spectral density of x(t) 

frequency, radians/second 
incremental value of c o, radians/second 
estimate 


Matrix notation: 

(x), x column matrix 

[X] rectangular or square matrix 

XT transpose 

X”1 inverse 

* complex conjugate 

Numbers used as subscripts denote the pertinent term or terms of the Volterra integral 
series or summation. 
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Figure 13. — Effect of smoothing on cross-spectra estimates 
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